dblogr/

Vaccine Hesitancy in Canada

Graphs of vaccine hesitancy in Canada using STATCAN data


Prepare Data

# devtools::install_github("derekmichaelwright/agData")
library(agData)
library(readxl) # read_xlsx()
# Prep data
myCaption <- "www.dblogr.com/ or derekmichaelwright.github.io/dblogr/ | Data: STATCAN"
myColors <- c("darkgreen", "darkred", "steelblue")
#
d1 <- read_xlsx("data_canada_vaccines.xlsx", "Table 1") %>%
  mutate(Measurement = factor(Measurement, levels = unique(.$Measurement)),
         Type = factor(Type, levels = unique(.$Type)))
d2 <- read_xlsx("data_canada_vaccines.xlsx", "Table 2") %>%
  mutate(Measurement = factor(Measurement, levels = unique(.$Measurement))) %>%
  gather(Age, Value, 2:4)

By Group

# Plot
mp <- ggplot(d1, aes(x = Measurement, color = Color)) +
  geom_point(aes(y = Value)) +
  geom_errorbar(aes(ymin = Min, ymax = Max)) +
  facet_grid(. ~ Type, scales = "free_x", space = "free_x") +
  scale_color_manual(values = myColors) +
  scale_y_continuous(breaks = seq(50, 100, by = 10)) +
  theme(strip.background = element_blank(),
               strip.text.x = element_blank()) +
  theme_agData(legend.position = "none",
               axis.text.x = element_text(angle = 45, hjust = 1)) +
  labs(title = "Percent willing to receive the COVID-19 vaccine", 
       y = "Percent", x = NULL, caption = myCaption)
ggsave("canada_vaccines_01.png", mp, width = 6, height = 4)

By Age

# Plot
mp <- ggplot(d2, aes(x = Measurement, y = Value, fill = Age)) +
  geom_bar(stat = "identity", position = "dodge", 
           color = "black", alpha = 0.7) +
  scale_fill_manual(name = NULL, values = myColors[c(1,3,2)]) +
  theme_agData(legend.position = "bottom") +
  labs(title = "Percent willing to receive the COVID-19 vaccine", 
       y = "Percent", x = NULL, caption = myCaption)
ggsave("canada_vaccines_02.png", mp, width = 6, height = 4)


dblogr/


© Derek Michael Wright